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The form factors of weak decays of B s mesons to ground state D s mesons as 
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well as to their radial D s (2i>) and orbital D j excitations are calculated in the 
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framework of the relativistic quark model based on the quasipotential approach. 
£N| \ All relativistic effects, including contributions of intermediate negative-energy states 

C\ \ and boosts of the meson wave functions, are consistently taken into account. As 

pP a result the form factors are determined in the whole kinematical range without 

additional phenomenological parametrizations and extrapolations. On this basis 
semileptonic decay branching fractions are calculated. Two-body nonleptonic B s 
decays are considered within the factorization approximation. The obtained results 
' agree well with available experimental data. 
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I. INTRODUCTION 

> ■ 

^ , In recent years significant experimental progress has been achieved in studying properties 

of B s mesons [lj. The Belle Collaboration considerably increased the number of observed 
B s mesons and their decays due to the data collected in e + e~ collisions at the T( 10860) 
resonance 0. On the other hand, B s mesons are copiously produced at Large Hadron Col- 
lider (LHC). First precise data on their properties are coming from the LHCb Collaboration 
[3|. Several weak decay modes of the B s meson were observed for the first time jlj. New 
data are expected in near future [Bj]. The study of weak B s decays is important for further 
improvement in the determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix ele- 
ments, for testing the prediction of the Standard Model and searching for possible deviations 
from these predictions, the so-called "new physics" . 

The dominant decay channel of the B s meson is into the D s meson plus anything [lj]. 
Therefore various important properties of excited D s mesons can be studied in the B s meson 
weak decays. In particular, they can shed light on the controversial _D* (2317) and -D s i(2460) 
mesons, whose nature still remains unclear in the literature. The abnormally light masses 
of these mesons put them below DK and D*K thresholds, thus making these states narrow 
since the only allowed decays violate isospin. The peculiar feature of these mesons is that 
they have masses almost equal to or even lower than the masses of their charmed counterparts 
Z>2(2400) and £>i(2420) [1]. If these mesons are indeed 1 3 P and IP X states, then all IP states 
of the D s meson are narrow, contrary to the D meson case. This narrowness considerably 
simplifies the experimental investigation of weak B s decays to orbitally excited D sJ mesons. 
Recently it was proposed [sj that study of B s — > D^*] transitions can clarify some puzzles in 
the corresponding semileptonic B decays. 

In this paper, we extend our investigations of weak B and B c decays 0, Hj] to studying 
exclusive weak semileptonic and nonleptonic decays of the B s to the ground state, radially 
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and orbitally excited D s mesons. For the calculations we use the same effective methods 
0,0] previously developed and successfully applied in the framework of the QCD-motivated 
relativistic quark model based on the quasipotential approach. The weak decay matrix el- 
ements are parametrized by the invariant form factors which are then expressed through 
the overlap integrals of the meson wave functions. The systematic account for relativistic 
effects, including the wave function transformations to the moving reference frame and con- 
tributions from the intermediate negative-energy states, allows one to reliably determine the 
momentum transfer dependence of the decay form factors in the whole accessible kinematical 
range. The other important advantage of this approach is that for numerical calculations we 
use the relativistic wave functions, obtained in the meson mass spectra calculations [9I. Ho|. 
Thus we do not need any additional ad hoc parametrizations or extrapolations which were 
usually used in some previous investigations. 

The calculated form factors are then substituted in expressions for the differential decay 
rates and semileptonic decay branching fractions are evaluated. The tree-dominated two- 
body nonleptonic B s decays to the D s meson and light or another D s meson are studied on 
the basis of the factorization approach. Such approximation significantly simplifies calcula- 
tions, since it allows one to express the matrix elements of the weak Hamiltonian governing 
the nonleptonic decays through the product of the transition matrix elements and meson 
weak decay constants. All these ingredients are available in our model. The obtained results 
are compared with previous calculations and experimental values, which are measured for 
some B s semi-exclusive semileptonic and several exclusive nonleptonic decay modes. 

The paper is organized as follows. In Sec. [IT] we briefly describe the relativistic quark 
model. Then in Sec. II III we discuss the relativistic calculation of the transition matrix el- 
ement of the weak b — > c current between meson states in the quasipotential approach. 
Special attention is paid to the contributions of the negative energy states and the relativis- 
tic transformation of the wave functions to the moving reference frame. These methods are 
applied in Sec. IIVI to the calculation of the form factors of weak B s decays to ground state 
D s mesons. The form factors are obtained as the overlap integrals of meson wave func- 
tions within the heavy quark expansion up to subleading order. It is shown that all heavy 
quark symmetry relations are explicitly satisfied. These form factors are used for evaluating 
semileptonic decay branching fractions in Sec. |V] The calculations of the form factors and 
semileptonic decay branching fractions for B s decays to radially excited D^(2S) mesons 

are presented in Sees. I VI I and IVHI within the same approach. In Sec. IVIIII the form factors 

(*) 

of weak B s decays to orbitally excited D sJ mesons are obtained. Semileptonic branching 
fractions for B s decays to orbitally excited D^*J mesons are given in Sec. HXl Finally, the 
two-body nonleptonic B s decays calculated within the factorization approximation are pre- 
sented in Sec. |2D All obtained results are confronted with previous calculations and available 
experimental data. Section |XI] contains the conclusions. The relations between the sets of 
weak form factors, the model independent HQET expressions for the form factors, and he- 
licity components of the hadronic tensor defined in terms of the form factors are presented 
in the Appendices. 
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II. RELATIVISTIC QUARK MODEL 

In the quasipotential approach a meson is described as a bound quark- ant iquark state 
with a wave function satisfying the quasipotential equation of the Schrodinger type |9|] 

where the relativistic reduced mass is 

E X E 2 M 4 - (mj - m 2 ) 2 
^ " E x + E 2 ~ 4M3 ' [Z) 

and Ei, E 2 are the center of mass energies on mass shell given by 

M 2 — m\+m\ M 2 — mf + m\ 

1_ 2M ' 2 ~ 2M ' ' ^ 

Here M = Ei + E 2 is the meson mass, mi^ are the quark masses, and p is their relative 
momentum. In the center of mass system the relative momentum squared on mass shell 
rGcicls 

h 2 (M) _ l M " ~ ( m i + m 2 ) 2 }[M 2 - (mi - m 2 ) 2 ] 
1 1 AM 2 ' 1 ' 

The kernel V(p, q; M) in Eq. ([T]) is the quasipotential operator of the quark-antiquark in- 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in- 
teraction, we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by [s| 

V(p, q; M) = Ui(p)u 2 (-p)V(p, q; M)ui(q)u 2 (-q), (5) 

with 

V(p, q; M) = ^a s D, u (k)^r 2 + V c v oni (k)T^ 2 ,, + V* ni (k), 
where a s is the QCD coupling constant, D^ v is the gluon propagator in the Coulomb gauge 

D°°(k) = -^, U«(k) = -^^-^V D* = D* = 0, (6) 

and k = p — q. Here 7^ and u(p) are the Dirac matrices and spinors 

u\v) 



t ^ ± m ( crp ) S (7) 

where a and x X are Pauli matrices and spinors, respectively, and e(p) = a/p 2 + m 2 . The 
effective long-range vector vertex is given by 

r M (k) = 7M + ^^*", ( 8 ) 
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where k is the Pauli interaction constant characterizing the long-range anomalous chromo- 
magnetic moment of quarks. Vector and scalar confining potentials in the nonrelativistic 
limit reduce to 

KonfW = (l-e)(Ar + B), 

Ki f (r) = e(Ar + B), (9) 

reproducing 

Konf (r) = Ko» + Konf (r) = Ar + B, (10) 

where e is the mixing coefficient. 

The expression for the quasipotential of the heavy quarkonia within and without the 
v 2 /c 2 expansion can be found in Ref. [9|. The quasipotential for the heavy quark interaction 
with a light antiquark without employing the nonrelativistic (v/c) expansion is given in 
Ref. All the parameters of our model like quark masses, parameters of the linear 

confining potential A and B, mixing coefficient e and anomalous chromomagnetic quark 
moment k are fixed from the analysis of heavy quarkonium masses and radiative decays 0] . 
The quark masses m h = 4.88 GeV, m c = 1.55 GeV, m s = 0.5 GeV, m U)d = 0.33 GeV and the 
parameters of the linear potential A = 0.18 GeV 2 and B = —0.30 GeV have values inherent 
for quark models. The value of the mixing coefficient of vector and scalar confining potentials 
e = — 1 has been determined from the consideration of the heavy quark expansion for the 
semileptonic B — > D decays [ll[ and charmonium radiative decays • Finally, the universal 
Pauli interaction constant k = — 1 has been fixed from the analysis of the fine splitting of 
heavy quarkonia 3 Pj- states 0] and the heavy quark expansion for semileptonic decays of 
heavy mesons [ll[ and baryons [13]. Note that the long-range magnetic contribution to the 
potential in our model is proportional to (1 + k) and thus vanishes for the chosen value of 
k = — 1 in accordance with the flux tube model. 



III. MATRIX ELEMENTS OF THE ELECTROWEAK CURRENT FOR b -> c 

TRANSITION 

In order to calculate the exclusive semileptonic decay rate of the B s meson, it is necessary 
to determine the corresponding matrix element of the weak current between meson states. 
In the quasipotential approach, the matrix element of the weak current Jjf = C7 M (1 — 75)6, 
associated with the b — > c transition, between a B s meson with mass Mg and momentum 



Pb s and a final D s meson with mass Md s and momentum pr> 3 takes the form [13 



(D.(pd.)\J?\B.(pb.)) = J ^^^P, s (p)r,(p,q)^ Pfls (q), (11) 

where r^(p, q) is the two-particle vertex function and ^m Pm are the meson (M = B S ,D S ) 
wave functions projected onto the positive energy states of quarks and boosted to the moving 
reference frame with momentum pm- 

The contributions to r^(p,q) come from Figs. [1] and [2j The contribution T^(p, q) is 
the consequence of the projection onto the positive-energy states. Note that the form of the 
relativistic corrections emerging from the vertex function rj, 2 ' (p, q) explicitly depends on the 
Lorentz structure of the quark-antiquark interaction. In the heavy quark limit mq — > 00, 



B s 



Pc 



FIG. 1: Lowest order vertex function rW(p, q) contributing to the current matrix element (jlip . 




<is . Ps <is . Ps 

FIG. 2: Vertex function r( 2 )(p,q) taking the quark interaction into account. Dashed lines 
correspond to the effective potential V in ([5]). Bold lines denote the negative-energy part of the 
quark propagator. 



only r^(p,q) contributes, while T^(p, q) give contributions starting from the subleading 
order. The vertex functions look like 



r«(p, q) = u c (p c ) 7M (l - 7 5 H(g 6 )(27r) 3 5(p, - q s 



(12) 



and 



r i 2) (p ; q) = u c (p c )u s (p s ){V(ps - qs 



+7i,(l -7r 5 ) Aif){k) ~° 



£b(k) + e b (p c ) 



e c {k') + e c {q b ) 

7?V(p s - - q s )}u b (q b )u s (q s ) } 



(13) 



where the superscripts "(1)" and "(2)" correspond to Figs. [T] and [2} k = p c — A; k' 
q& + A; A = p Ds - p Bs ] 



AH(p) 



e(p) - (m7° + 7 (7p)) 
2^ 



Here 13 



PCS = ZcM^±Y,^\vD a )v\ 
M D S i= l 

M B a i=l 
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and are three four-vectors given by 

n®»(p) = 1^7, 6a + „, p y 7tv 1 , E = \/p 2 + M 2 . 
\M' j M(£ + M)J' vh 

The wave function of a final Z) s meson at rest is given by 

*D S ( P ) = = 3^ ^ sJ ( P ), (14) 

where J and .M are the total meson angular momentum and its projection, L is the orbital 
momentum, while S — 0, 1 is the total spin. ipn sJ (p) is the radial part of the wave function, 
which has been determined by the numerical solution of Eq. (JTJ) in Ref. [l(|. The spin-angular 
momentum part y JLSM has the following form 

y jLSM = J2(LM-ax-a 2 , S a x +a 2 \J M){\a x , \a 2 \S a x + ^Yf-^xiMxifa). 

<Tl<T 2 

(15) 

Here (71J711, j2 m 2|L.M) are the Clebsch-Gordan coefficients, YJ m are the spherical harmon- 
ics, and x{°) (where a = ±1/2) are the spin wave functions, 

X(V2)= Q, X(-1/2) = (J). 

The heavy-light meson states (D s i, D' sl ) with J — L — 1 are the mixtures of spin-triplet 
_D s ( 3 Pi) and spin-singlet D s ( 1 Pi) states: 

= ^^(ipocos^ + ^D^sp^sinv?, 
ty D > si = -i&D^Pt) sin y> + i&D.pPi) cos (16) 

where = 34.5° is the mixing angle and the primed state has the heavier mass [Toj]. Such 
mixing occurs due to the nondiagonal spin-orbit and tensor terms in the Qq quasipotential. 
The physical states are obtained by diagonalizing the corresponding mixing terms. Note 
that the above value of the mixing angle ip is very close to its heavy quark limit ip mQ ^oo = 



arctan(yl/2) rs 35.3°. This means that the wave functions ^d s1 and correspond in 
the heavy quark limit to ^d s(3/2) and ^ d s{1/2) , respectively. 

It is important to point out that the wave functions entering the weak current matrix 
element (TTTj) are not in the rest frame in general. For example, in the B s meson rest frame 
(Pb s = 0), the final meson is moving with the recoil momentum A. The wave function of 
the moving meson o s a i s connected with the wave function in the rest frame ^ d s o = ^ d s 



by the transformation [13 



* D .a(p) = ^WJ^WJWp), (17) 

where R w is the Wigner rotation, La is the Lorentz boost from the meson rest frame to a 
moving one, and the rotation matrix D l l 2 (R) in spinor representation is given by 

(J fjD 1 s ! c 2 (RYj = S-\p s , c )S(A)S(p) } (18) 

where 



s(p) = \ e{p l + m [1+ ap 



2m \ e{p) + m J 
is the usual Lorentz transformation matrix of the Dirac spinor. 
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IV. FORM FACTORS OF WEAK B s DECAYS TO D s MESONS 

For considering weak B s decays to ground state D s mesons we employ the heavy quark 
expansion which significantly simplifies calculations. Therefore it is convenient to intro- 
duce the heavy quark effective theory (HQET) parametrization for the weak decay matrix 
elements 14, 15| : 



{D.V)\*fW.V)) h+{v + vT + h _ {v ^ vr , (19) 

0, (20) 
ihye^elvfa, (21) 

h Al (w + - (h M v» + h A3 v'»)(e* ■ v), (22) 

where v (?/) is the four- velocity of the B s (Dy)) meson, e M is the polarization vector of the 
final vector meson, and the form factors hi are dimensionless functions of the product of 
four-velocities 

W = V ■ V = 





^M Ds M Ba 


(D s ( 


v')\^b l5 \B s {v)) 


(D* s ( 


v',e)\ry%\B s (v)) 




\J Md*Mb 3 




,e)\cri 5 b\B s (v)} 


jM Dt M Bs 



2M Bs M r 



,(*) 



and q = Pb s — P n (*) is the momentum transfer from the parent to daughter meson, M Bs is 
the B s meson mass, M „(*) is the final meson mass and e„ is the polarization vector of 
the final vector D* meson. 

In HQET these form factors up to 1/ mg order are expressed through one leading Isgur- 



Wise function £, four sudleading functions Xi,2,3 and one mass parameter A lq . These 
relations are given in Appendix [A] 

To calculate the weak decay matrix element in the quasipotential approach, we substitute 
the vertex functions (fT2|) and (fl3l) in Eq. (TlT!) and take into account the wave function 
transformations (fT7|) . The contribution of the leading order vertex function r^(p, q) can 
be easily simplified by carrying out one of the integrations using the 5-function. Then we 
employ the heavy quark expansion, which permits us to take one of the integrals in the 
contribution of the vertex function T^(p, q) to the weak current matrix element. As a 
result we express all matrix elements through the usual overlap integrals of the meson wave 
functions. We carry out the heavy quark expansion up to the second order and compare the 
obtained expressions with model independent HQET relations flAlj) - (IA6p . 

All leading order relations are exactly satisfied. In this limit of an infinitely heavy quark, 
all form factors are expressed through the single universal Isgur-Wise function £(w) |l4| 

h + (w) = h Al (w) = h Aa (w) = h v (w) = £(w) 

h-(w) = h M {w) = 0. (23) 



This function is given by the following overlap integral of meson wave functions [11 



£H = \hrrrJ^ I ^*D s [p + 2e s (p)J——e± I * fls (p), (2 I) 



W + 1 m O- > °° J (2lT 
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where = A/V A 2 is the unit vector in the direction of A = M^v' — Me a v. In the 
infinitely heavy quark mass limit the wave functions of initial tyft s and final ^£> s heavy 
mesons coincide. As a result the HQET normalization condition [15( 

£(1) = 1 

is exactly reproduced. 

In order to fulfill the HQET relations (jAip -( TA~6]) at the first order of the heavy quark 
I/tuq expansion it is necessary to set (1 — e)(l + k) = 0, which leads to the vanishing long- 
range chromomagnetic interaction. This condition is satisfied by our choice of the anomalous 
chromomagnetic quark moment k — — 1. To reproduce the HQET relations at second order 
in l/m,Q, one needs to set e = —1 [llj. This serves as an additional justification, based 
on the heavy quark symmetry and heavy quark expansion in QCD, for the choice of the 
characteristic parameters in our model. The subleading Isgur-Wise functions are given by 

Uw) = (A — m q ) (l + £H, (25) 

XiH = A^l±£H, (26) 
w + 1 

= -S^TT« ra) ' (27) 

X3{w) = Te^^Tl^' (28) 

where the HQET parameter A = M — itlq is equal to the mean energy of a light s quark in 
a heavy meson (e s ). The functions Xi an d X3 explicitly satisfy normalization conditions at 



the zero recoil point 16 



Xi(l) = Xs(l) = 0, 

arising from the vector current conservation. 

Near the zero recoil point of the final meson w — 1 the Isgur-Wise functions have the 
following expansions 

£{w) = 1 - 1.466(w - 1) + 1.844(w - l) 2 + • • • , 

£a(w)/A = 0.359 - 0.408(u; - 1) + 0.428(u> - l) 2 + ■ ■ ■ , 

XiH/A = 0.499(w - 1) - 0.982(u; - l) 2 + ■ ■ • , 

X2H/A = -0.0156 + 0.0307(w - 1) - 0.0442(w - l) 2 + • • ■ , 

XsH/A = 0.0312(w - 1) - 0.0614(w - l) 2 + ■ • • (29) 

The calculated leading and subleading Isgur-Wise functions for B s — > transitions 
are plotted in Fig. Using these Isgur-Wise functions we obtain the decay form factors 
hi(w) with the account of the first order l/mq corrections in the whole kinematical range. 
To improve calculations we scale the results by the values of the form factors at zero recoil 



hi(l) evaluated with the inclusion of 1/toq corrections using formulas given in Ref. [llj. We 
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FIG. 3: Leading and subleading Isgur-Wise functions for B s — )■ D s transitions. 



find that the account of the first order corrections changes the form factors by about 17%, 
while the contribution of the second order corrections is less than 3%. These values are in 
accord with the naive estimates of such corrections A/(2m c ) ~ 0.15 and [A/(2m c )] 2 0.02. 

The other popular parametrization for the matrix elements of weak current J w between 
meson states is given by 



(D s (p Ds )\c^b\B s (p Bs )) = f + {q 2 ) 



M% -M 2 D ' 
Pb s + Pd s 2 q 



„ , 9 , Ml - M 2 D , 
+ fo(q 2 ) Bs „ 9 Ds 



r 



(D* s (p D *)\c^b\B(p Bs )) 



(D s ( PDs )\c^ l5 b\B s ( PBs )) = 0, 

2iV(q 2 ) 
M Bs + M D . 



e^ pa e; PBapPD:cT , 



(30) 
(31) 



(32) 



(Dt( PD *)\crjsb\B s ( PBs )) = 2M DS M^' ' " 



q^ + (M Bs + M^JA^q 2 ) [e*» - 



e ■ q 



~A 2 (q 2 



e ■ q 



PBs + Pot ~ 



Ml - M 2 



(33) 



M Ba + M D% 

At the maximum recoil point (q 2 = 0) these form factors satisfy the following conditions: 

/+(0) = / (0), 
, . , M B + M D * , . . M B — M B * . , s 

The relations between two sets of weak decay form factors are given in Appendix Q31 
Substituting in these relations the Isgur-Wise functions of our model (124"]) - (128]) we find 

that the decay form factors can be approximated with sufficient accuracy by the following 

expressions: 

(a) f + (q 2 ),V(q 2 ),A (q 2 ) = F(q 2 ) 



F{q 2 



F(0) 



M 2 



' q 2 q 4 



(34) 



M 



b: 



B: 
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TABLE I: Form factors of weak B s — >■ D s transitions calculated in our model. Form factors 
f+(Q 2 ), V(q 2 ), A (q 2 ) are fitted by Eq. ([Ml), and form factors f (q 2 ), A^q 2 ), A 2 {q 2 ) are fitted by 



Eq. {35]). 
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FIG. 4: Form factors of weak — )• D s transitions. 



(b) f (q 2 ),A 1 (q 2 ),A 2 (q 2 )=F(q 2 ) 

Hi 2 ) = t ^ rr» ( 35 ) 

where M = M B; = 6.332 GeV for the form factors f+(q 2 ), V(q 2 ) and M = M Bc = 6.272 GeV 
for the form factor Ao(q 2 ); the values F(0) and a\ )2 are given in Tabled! The values of a\^ 2 are 
determined with a few tenths of percent errors. The main uncertainties of the form factors 
originate from the account of 1/mg corrections at zero recoil only and from the higher order 
I /trig contributions and can be roughly estimated in our approach to be about 2%. 1 The 
q 2 dependence of these form factors is shown in Fig. HI 



Other uncertainties originating, e.g., from meson wave functions and model parameters are significantly 
smaller. Indeed meson wave functions and masses were obtained by the numerical solution of the quasipo- 
tential equation with the completely relativistic spin-independent and spin-dependent potentials treated 
nonperturbatively [lo| . The model parameters were fixed in previous calculations which correctly repro- 
duce numerous experimental data. The integrated quantities such as decay form factors and semileptonic 
decay rates are much less sensitive to the variation of the model parameters than such quantities as hadron 
masses which are measured with considerably higher accuracy. Thus even the limited variation of these 
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TABLE II: Comparison of theoretical predictions for the form factors of semileptonic decays B s — > 
Ds*^eu at maximum recoil point q 2 = 0. 

Mo) vxo) Mo) Mo) a 2 (o) 

this paper 0.74 ±0.02 0.95 ± 0.02 0.67 ±0.01 0.70 ± 0.01 0.75 ± 0.02 
[17] 0.61 0.64 0.56 0.59 

[18] 0.7 ±0.1 0.63 ±0.05 0.52 ± 0.06 0.62 ± 0.01 0.75 ± 0.07 

[19] 0.57iH! 0.70iHI 0.651H1 0.6^ 

[20] 0.86±g:g 



[21] 0.74i°;°5 o.63l8;81 0.61±8;°t 0.59". 



-0.04 
•0.04 



In Table HI] we confront our predictions for the form factors of semileptonic decays B s — > 
D^ev at maximum recoil point q 2 = with results of other approaches [l7l-21]. Different 
quark models are used in Refs. [l7|, 19, 2l||, while the QCD and light cone sum rules are 
employed in Refs. 18|, |20j . We find that these significantly different theoretical calculations 
lead to rather close values of the decay form factors. One of the main advantages of our 
model is its ability not only to obtain the decay form factors at the single kinematical 
point, but also to determine its q 2 dependence in the whole range without any additional 
assumptions or extrapolations. 



V. SEMILEPTONIC B s DECAYS TO D s MESONS 

The differential decay rate for the semileptonic B s meson decay to mesons reads 22| 



dT(B s ^D^W) G 2 F 2 A 1/2 (g 2 - mf) 2 

/ d>\ 



(2tt 



24M| q 2 



(36) 



where Gf is the Fermi constant, V c b is the CKM matrix element, A = A(M| a , M^*), q 
M% a + Af* w + <? 4 - 2{M%M 2 dW + M^q 2 + M| s g 2 ), m ; is the lepton mass and 



HH ] = H + Hl + H_Hl_ + H Hl 



(37) 



The helicity components H±, H and H t of the hadronic tensor are expressed through the 
invariant form factors. They are given in Appendix O 

Now we substitute the weak decay form factors calculated in the previous section into the 
above expressions for decay rates. The resulting differential decay rates for the B s decays 
to the mesons are plotted in Fig. [51 The corresponding total decay rates are obtained 
by integrating the differential decay rates over q 2 . For calculations we use the CKM matrix 
element | V c b\ = (3.9±0.15) x 10~ 2 , which was obtained from the comparison of our theoretical 
predictions 0, 23] for the products F D (*)(w)\V c b\ and for the B — > D^Hui decay branching 
fractions with updated experimental data. 2 It is necessary to point out that the kinematical 



parameters, permitted by the description of hadron masses, will give significantly smaller contributions 
to the form factor and decay rate uncertainties compared to the ones mentioned above. 
2 This value of \V c b\ is in accord with its recent evaluation by the Heavy Flavor Averaging Group [24j . 
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FIG. 5: Predictions for the differential decay rates of the B s — >■ D s Iv semileptonic decays. 



TABLE III: Comparison of theoretical predictions for the branching fractions of semileptonic decays 
B s D^lv (in %). 





Decay 


this paper 


[18] 


[19] 


[20] 


[21] 


[25] 


[26] 


B & 


—> D s ev 


2.1 ±0.2 


1.35 ±0.21 


1.4-1.7 






2.73-3.00 


2.8-3.8 


B s 


-> Z? s rz^ 


0.62 ±0.05 




0.47-0.55 


o.33±8;i! 








B s 


-»■ D* s ev 


5.3 ±0.5 


2.5 ±0.1 


5.1-5.8 




5.2 ±0.6 


7.49-7.66 


1.89-6.61 


B s 




1.3 ±0.1 




1.2-1.3 




i q+0.2 







range accessible in these semileptonic decays is rather broad. Therefore the knowledge of the 
q 2 dependence of the form factors is very important for reducing theoretical uncertainties 
of the decay rates. Our results for the semileptonic B s — > DyHu decay rates are given in 
Table II III in comparison with previous calculations. The authors of Ref. 18| use the QCD 



sum rules, while the light cone sum rules approach is adopted in Ref. [20J. Different types 



of constituent quark models are employed in Refs. [19|, |21|, |25| and the three point QCD 
sum rules are used in Ref. (26|. We see that our predictions are consistent with results of 



quark model calculations in Refs. [19|, |21] . They are approximately two times larger than 
the QCD sum rules an d lig ht cone sum rules results of Refs. [H, 2(| , but slightly lower than 
the values of Refs. 25, 26 



We find that the total branching fraction of the semileptonic decays of B s mesons to the 
ground state £>W i s equal to Br(B s -»■ £>|*W) = (7.4 ± 0.7)% and Br(B s -»■ D^ru) = 
(1.92 ± 0.15)%. The errors in our estimates originate from the uncertainties in the deter- 
mination of the CKM matrix element \V c b\, which are dominant, and from the theoretical 
uncertainties in the determination of decay form factors. The latter uncertainties are consid- 
erably smaller than the former ones and are mostly related with the estimates of the higher 
order terms in the heavy quark expansion. 



VI. FORM FACTORS OF WEAK B s DECAYS TO RADIALLY EXCITED D S (2S) 

MESONS 

The decay form factors (IT9"l) - (j2"2"l) up to 1/mq order in HQET for B s decays to radially 
excited D s [(n ±1)5"] mesons are expressed through one leading and five sub leading £ 3 , 
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Xi,2,3,f> Isgur-Wise functions and two mass parameters A and A^ n ^ 22|. They are presented 
in Appendix IDl 



In our model all HQET relations ( 1Dll) -( jD6l) are satisfied and we get the following expres 



sions for the leading and subleading Isgur-Wise functions [2 



| 3 „ = (^±1 _ m , + l^A) U + (m(w) , (39 ) 



6 w — 1 J \ 3w + 1 

1 u> - 1 AW - A m . 
£ w (w) 

^ / (^^(p + 2 ,W^e a ) ,<?>, P ), ,40) 

S -l2^1^ f( "W- (41 > 
80 it; + 1 it; — 1 



where A 2 = M 2 W (w 2 — 1). Here we used the expansion for the S'-wave meson wave 

r 1 D s (2o) 

function [27( 

^ M = Vff + Aa/^q (yfi* + d M ^ sd ) + Oil/ml), 

where is the wave function in the limit uiq — > oo, ari d ^ are the s P m_ 

independent and spin-dependent first order 1/rriQ corrections, dp = — 3 for pseudoscalar 
and dy — 1 for vector mesons. The symbol = in the expressions fj4*0"]) - (j4*3l for the sublead- 
ing functions Xi(w) implies that corrections suppressed by an additional power of the ratio 
(w — l)/(w + 1), which is equal to zero at w — 1 and less than 1/6 at w max , were neglected. 
Since the main contribution to the decay rate comes from the values of form factors close 
to w — 1, these corrections turn out to be unimportant. 

It is clear from the expression (|38|) that the leading order contribution vanishes at the 
point of zero recoil (A = 0, w = 1) of the final D^(2S) meson, since the radial parts of the 
wave functions *&d s (2S) an d ^b s are orthogonal in the infinitely heavy quark limit. 

Near the zero recoil point of the final meson w = 1 the Isgur-Wise functions have the 
following expansions 

£ (1) (u0 = 2.455(w - 1) - 9.545(w - l) 2 + ■ ■ • , 
fsW/A = 0.221 + 1.094O - 1) - 5.294(iu - l) 2 + ■ ■ • , 
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FIG. 6: Leading and subleading Isgur-Wise functions for B s — > (25) transitions (A = (A^ + 
A)/2). 



TABLE IV: Form factors of weak B s — > (25) transitions calculated in our model 





B s - 


-> D S (2S) 




B s - 


*D* S (2S) 






f+ 


fo 


V 


A 


At 


A 2 


F(0) 


0.41 


0.41 


0.46 


0.45 


0.32 


0.058 


-^(?max) 


0.28 


0.075 


0.28 


0.29 


0.072 


-0.29 



XiH/A = 0.182 - 0.143(w - 1) - 1.055O - l) 2 + ■ ■ • , 

X2H/A = -0.0552 + 0.242(w - 1) - 0.547(w - l) 2 + ■ ■ • , 

XsH/A = -0.00133 - 0.0334(u; - 1) + 0.150(w - l) 2 + • • ■ , 

X b (w)/A = -0.169 + 1.114(w - 1) - 4.060(w - l) 2 + ■ ■ • , (44) 

where A = (A (1) + A)/2. 

The leading and subleading £ 3 , Xi,2,3,b Isgur-Wise functions for B s — > D^(2S) tran- 
sitions are shown in Fig. [6j We use relations (1Blj) - (1B6j) to express form factors f +i0 (q 2 ), 
V(q 2 ) and Ad, 1,2 through the calculated Isgur-Wise functions. The obtained form factors 
are plotted in Fig. [71 Their values at zero and maximum q 2 are given in Table IIVI The 
main theoretical uncertainties of the decay from factors, as for the decays to the ground 
state mesons, originate from the higher order 1/mg contributions and are less than 4%. 
Comparing plots in Figs. H] and [7J we see that form factors for the decays to ground and 
radially excited states have significantly different behaviour in q 2 . The former ones grow 
with q 2 , while the latter ones decrease. This is the consequence of the different structure of 
nodes of the wave functions of these states. 
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FIG. 7: Form factors of the weak B s — > Di (2S) transitions. 



TABLE V: Predictions for the branching fractions of semileptonic decays B s — > Dg (2S)lv (in %) 





Decay 


Br 


B s 


-»■ D s {2S)ev 


0.27 ±0.03 


B s - 


D s (2S)tu 


0.011 ±0.001 


B s - 


-> D*(2S)eu 


0.38 ± 0.04 


B s - 


^D* s {2S)tv 


0.015 ±0.002 



VII. SEMILEPTONIC B s DECAYS TO RADIALLY EXCITED D S {2S) MESONS 



For the calculation of the semileptonic B s decays to radially excited D S (2S) mesons 
we use the expression for the differential decay rates ( 136|) with the helicity components 
of the hadronic tensor given by Eqs. ( lCip -( IU4|) and decay form factors calculated in the 
previous section. The predictions for the corresponding branching fractions are given in 
Table [V] We find that semileptonic B s decays to the pseudoscalar D S (2S) and vector D*(2S) 
mesons have close values. The total contribution of these decays is obtained to be Br(B s — > 
D^;\2S)ev) = (0.65 ± 0.06)% and Br{B s -> D^{2S)tv) = (0.026 ± 0.003)%. 




01234567 01234567 



q 2 (GeV 2 ) q 2 (GeV 2 ) 



FIG. 8: Predictions for the differential decay rates of the B s — > D s (2S)lu semileptonic decays. 
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The differential decay rates of the B s — > {2S)lv semileptonic decays are plotted in 
FigEl 

VIII. FORM FACTORS OF WEAK B s DECAYS TO ORBITALLY EXCITED D { *] 

MESONS 

The matrix elements of the weak current Jj^ = ^7^(1 — 7s)c for B s decays to orbitally 
excited scalar D* mesons can be parametrized by two invariant form factors 

(D* s0 (p Ds0 )\crb\B s (p Bs )} = 0, (45) 
(D: (p Ds0 )\cYlsb\B(p Ba )) = r + (q 2 ) (p Ba + pQ +r_(g 2 ) (p Ba ~ pQ , (46) 

where q = p Ba — Pd s0 , ^d s0 is the scalar meson mass. 

The matrix elements of the weak current for B s decays to the axial vector D s i meson can 
be expressed in terms of four invariant form factors 

(D sl (p Dal )\c^b\B(p Ba )) = (M Ba + M Dsl )h Vl (q 2 )e*» + [hy 2 {q 2 )p Ba + M^KJ^r 1 , (47) 

M B S 

(D s i(PD 3l )\cj^b\B(p Bs )) = J^fff e^e: PBsP p Dsl<7 , (48) 

where Md s1 and e M are the mass and polarization vector of the axial vector meson. The 
matrix elements of the weak current for B s decays to the axial vector D' sl meson are obtained 
from Eqs. (14TI) by the replacement of the set of form factors hi(q 2 ) by gi(q 2 ) (i = Vi, V 2 , V 3 , A). 

The matrix elements of the weak current for B s decays to the tensor D* 2 meson can be 
decomposed in four Lorentz-invariant structures 

(D; 2 (p Ds2 )\crb\B(p Ba )} = J ttv{q J e^e: a §^ PBapPDs2a , (49) 

M Ba + M Da2 M Ba 

(D: 2 (p Da2 )\qr l5 b\B(p Ba )) = (M Ba + M Da2 )t Al (q 2 )e*^^ 

+ [tA 2 (q 2 )p Ba + tA 3 (q 2 )p% a2 Kp?j^, ( 5 °) 

where Md s2 and are the mass and polarization tensor of the tensor meson. 

(*) 

To obtain the form factors of B s — > D sJ weak transitions we use the expression for the 



weak current matrix element (TTTj) . We calculate exactly the contribution of the leading 
vertex function T^^p, q) ( [12]) to the transition matrix element of the weak current ( TTTj) 
using the 5-function. For the evaluation of the subleading contribution r^ 2 ^(p, q) for the 
B — > D^*J transitions, governed by the heavy-to-heavy b — > c transitions, we use expansions 

in inverse powers of masses of the heavy b- and c-quarks, contained in the initial B s meson 

(*) 

and final D sJ meson. Thus we can neglect the small relative quark momentum |p| compared 
to the heavy quark mass ttiq in the quark energy e<g(p + A) = Jtoq + (p + A) 2 , replacing it 



by cq(A) = Jtoq + A 2 in expressions for the r^ 2 ^(p, q). Note that we keep the dependence 
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FIG. 9: Form factors of the B s 


decays to the P-wave D s j mesons. 







on the recoil momentum A = p (,) — p_e s . This replacement removes the relative momentum 

dependence in the quark energy and thus permits us to perform one of the integrations in the 
r^ 2 -* (p, q) contribution using the quasipotential equation. The subleading contribution turns 
out to be rather small numerically, since it is proportional to the quark binding energy in the 
meson. Therefore we obtain reliable expressions for the form factors in the whole accessible 
kinematical range. It is important to emphasize that when doing these calculations we 
consistently take into account all relativistic contributions including boosts of the meson 
wave functions from the rest reference frame to the moving ones, given by Eq. ( TT7|) . The 
obtained expressions for the decay form factors are rather cumbersome and are given in the 
Appendix of Ref. [q. Note that, while calculating form factors of weak B s decays to D s \ 
and D' sl mesons, it is important to take into account the mixing (JIB!) of singlet D s QPx) and 
triplet D s ( 3 Pi) states. 

In Fig. M we plot form factors of the weak B s transitions to the P-wave D s j mesons. 
The calculated values of these form factors at q 2 = and q 2 = g^ ax = (M Bs — M DaJ ) 2 
are displayed in Table IVI1 The theoretical uncertainties of these form factors within our 
approach are mainly determined by the errors introduced by the replacement of e c (p + A) 
by e c (A) in the subleading vertex r^(p, q) and 0(l/mf) contributions. They are almost 
negligible at q 2 = and are less than 1% at q 2 = g^ax- 
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TABLE VI: Calculated values of the form factors of the B s decays to the P-wave D [ *j at q 2 = 
and q 2 = q 2 mSLX =(M Bs -M DsJ ) 2 . 





B s - 


^ u sG 


B s 








P. s ^ 


An 






5, -f A*2 






r+ 


r_ 


9A 3Vi 


9V 2 


9v 3 


h A 








tv 









0.40 


-0.91 


-0.61 0.01 


-0.24 


0.77 


-1.14 


-0.19 


-0.43 


-0.29 


-0.63 


-0.59 0.16 


-0.01 


q 2 


0.50 


-1.58 


-0.97 0.03 


-0.04 


1.29 


-1.77 


0.03 


-0.81 


-0.48 


-1.09 


-0.72 0.12 


-0.02 



TABLE VII: Comparison of the predictions for the branching fractions of the semileptonic decays 
B s Djjli/ (in %). 





Decay 


this paper 


m — )• oo 
[28] 


with 1/rriQ 
[28] 


[29] [25] 


[2Q] 


[30J 


m 


B s 


-> D* ez/ 


0.36 ±0.04 


0.10 


0.37 


0.443 0.49-0.571 


0-2318:13 


~ 0.1 


0.20 


B s 




0.019 ±0.002 










~ 0.01 




B s 




0.19 ±0.02 


0.13 


0.18 


0.174-0.570 0.752-0.869 




0.49 


0.10 


Bs 


-> ^ii^ 


0.015 ±0.002 














B s 


—7- D s iev 


0.84 ±0.09 


0.36 


1.06 


0.477 








B s 


-> AslTI/ 


0.049 ± 0.005 














B s 




0.67 ±0.07 


0.56 


0.75 


0.376 








B s 


-> A*2™ 


0.029 ± 0.003 















IX. SEMILEPTONIC B s DECAYS TO ORBITALLY EXCITED D { *] MESONS 

The differential semileptonic decay rates of B s mesons to orbitally excited d[*J mesons 
are given by Eq. fl36l) . The helicity components H±, H and H t of the hadronic tensor 
are expressed through the invariant form factors (j43|) - fl50|) by the relations [s| given in 
Appendix [El 

Substituting calculated form factors in these expressions we get predictions for the branch- 
ing fractions of the semileptonic B s decays to orbitally excited D s mesons. We find that 
decays to D s ± and D* 2 mesons are dominant. The obtained results are given in Table IVHI 
in comparison with other calculations. First we compare with our previous calculation [28 



which was performed in the framework of the heavy quark expansion. We give results found 
in the infinitely heavy quark limit (tuq — > oo) and with the account of first order I/tuq 
corrections. It was argued 0, [32| that I/tuq corrections are large and their inclusion sig- 
nificantly influences the decays rates. The large effect of subleading heavy quark corrections 
was found to be a consequence of the vanishing of the leading order contributions to the 
decay matrix elements, due to heavy quark spin-flavour symmetry, at the point of zero re- 
coil of the final charmed meson, while the subleading order contributions do not vanish at 
this kinematical point. Here we calculated the decay rates without application of the heavy 
quark expansion. We find that nonperturbative results agree well with the ones obtained 
with the account of the leading order 1/mq corrections [28(|. This means that the higher 
order in 1/mq corrections are small, as was expected. Then we compare our predictions 



with the results of calculations in other approaches. The authors of Refs. [25|, |29| employ 
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FIG. 10: Predictions for the differential decay rates of the B — > D jlv semileptonic decays. 



B s ->D ,1 V; 




different types of constituent quark models for their calculations. Light cone and three 
point QCD sum rules are used in Refs. [20|, |30j , while HQET and sum rules are applied in 
Ref. [3 if . In general we find reasonable agreement between our predictions and results of 
Refs. [3, HE , but results of the quark model calculations jUJ are slightly larger. 



In Fig. [TU] we plot the differential decay rates of the B — > D^*)lv semileptonic decays. 
The total semileptonic decay branching fractions to orbitally excited D s mesons are found 
to be Br(B s -> D { *)ev) = (2.1 ± 0.2)% and Br(B s -> D { *Jtu) = (0.11 ± 0.01)%. 

The first experimental measurement of the semileptonic decay B s — > D s \[iv was done by 
the DO Collaboration {33|. The branching fraction was obtained by assuming that the D s \ 
production in semileptonic decay comes entirely from the B s decay and using a prediction 
for Br(D sl -> D*K° S ) = 0.25. Its value 

Br(B s -> D sl Xfxu) D0 = (1.03 ± 0.20 ± 0.17 ± 0.14)% 



is in good agreement with our prediction 0.84 ± 0.9 given in Table IVHI 

Recently the LHCb Collaboration [34| reported the first observation of the orbitally ex- 
cited D* 2 meson in the semileptonic B s decays. The decay to the D s i meson was also 
observed. The measured branching fractions relative to the total B s semileptonic rate are 

Br{B s -> D* s2 Xfiu)/Br(B s -> X F ) LHC b = (3.3 ± 1.0 ± 0.4)%, 

Br{B s -> D sl Xfiu)/Br(B s -> X/iz/) LHC b = (5.4 ± 1.2 ± 0.5)%. 
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The D* 2 /D s i event ratio is found to be 

Br(B s -> D* s2 Xfiv)/Br(B s -> D s iX/iz/) LHC b = 0.61 ± 0.14 ± 0.05. 

These values can be compared with our predictions if we assume that decays to D s i and 
D* 2 mesons give dominant contributions to the ratios. Summing up the semileptonic B s 
decay branching fractions to ground state, first radial and orbital excitations of D s mesons, 
presented in Sees. [V] IVHl HXf we get for the total B s semileptonic rate Br(B s — > X\iv) = 
(10.2 ± 1.0)%. Then using the calculated values from Table [VTT1 we get 

Br(B s -> D* s2 /jis)/Br(B s -> X/iz/) theor = (6.5 ± 1.2)%, 

Br{B s D slf iv)/Br(B s X/iz/) theor = (8.2 ± 1.6)%, 

and 

Br(B s ->• D* s2 nis)/Br(B s -> D slt iiy) theoT = 0.79 ± 0.14. 

The predicted central values are larger than experimental ones, but the results agree with 
experiment within 2a. 



X. NONLEPTONIC DECAYS 

In the standard model nonleptonic B s decays are described by the effective Hamiltonian, 
obtained by integrating out the Vy-boson and top quark. For the nonleptonic B s decay to 
the ground state or excited D s meson and light meson governed by b — > c transition the 
effective Hamiltonian is given by [35[ 

H cS = %V*V uq [cMO^ + c 2 (/i)0 2 u ] , (51) 

where q — d,s. For the nonleptonic B s decay to two charmed mesons the effective Hamilto- 
nian (AB = 1) [35j reads 

Q 10 

H« = ^av eq Y. c *fa)°i- ( 52 ) 



The Wilson coefficients are evaluated perturbatively at the W scale and then are 
evolved down to the renormalization scale fi ~ mj by the renormalization-group equations. 
Functions { are the local four-quark operators. The tree level operators have the form 

0\ = (jbc)v-A(<fQ)v-A, 

0\ = [bjCi) V -A(<fjQi)v-A, (53) 

while the functions Oi (i = 3, .., 10) are the penguin operators. The following notations are 
used 

{w)vtA = (ftnO-TlsW- 

The amplitude of the nonleptonic two-body B s decay to D s and light M mesons can be 
expressed through the matrix element of the effective weak Hamiltonian H e g in the following 

way 



G p 



M(B S D S M) = (D s M\H cS \B s ) = j V^V ug [ Cl (D,M\0?\B s ) + c 2 (D s M\0^\B s )] 



(54) 
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The factorization approach, which is widely used for the calculation of two-body nonlep- 
tonic decays, such as B s — » D S M, assumes that the nonleptonic decay amplitude reduces to 
the product of a meson transition matrix element and a weak decay constant [361 ] . Clearly, 
this assumption is not exact. However, it is expected that factorization can hold for ener- 
getic decays, where one final meson is heavy and the other meson is light and energetic (13]. 
A more general treatment of factorization is given in Ref. (38l |. 

Then the B s — > D~M + decay amplitude can be approximated by the product of one- 
particle matrix elements. The matrix element (q = d, s) is given by 

(M+DjfaO? + c 2 0%\B° s ) « a l {D:\{bc) v . A \B° s )(M + \{uq) v . A \Q) 1 (55) 

where the Wilson coefficients appear in the following linear combination 

a x = ci + -^-c 2 (56) 

and N r is the number of colors. For numerical calculations we use the values of Wilson 



coefficients given in Ref. [39 



The similar expression holds for B s — > decays [20j], namely 

10 / y*y \ 

(DjD;\J2c^)O<;\B s ) « U-^[a 4 + a 10 + r g (a 6 + a 8 )] 

i=l \ V cb V cs ) 

x (D; I (be) V „ A I B° s )(Dj\(cs) v _ A \0), (57) 

where the second term in brackets results from the penguin contributions, which are small 
numerically. The coefficients a 2n = c 2n + c 2n _i/iV c , and r q can be found, e.g., in Ref. [2p| . 

The matrix element of the weak current between vacuum and a final pseudoscalar 
(P) or vector (V) meson is parametrized by the decay constants fpy 

(P|gi7 M 75? 2 |0> = if P p^ (V\q lltl q 2 \Q) = e,M v f v . (58) 

The pseudoscalar fp and vector fv decay constants were calculated within our model in 
Ref. [40 . It was shown that the complete account of relativistic effects is necessary to 



get agreement with experiment for decay constants especially of light mesons. We use the 
following values of the decay constants: /„- = 0.131 GeV, f p = 0.208 GeV, fx = 0.156 GeV, 
f K * = 0.214 GeV, f Ds = 0.260 GeV and f D , = 0.315 GeV. The relevant CKM matrix 
elements are \V ud \ = 0.975, \V US \ = 0.225, \V CS \ = 0.973, \V ts \ = 0.0404, |y t6 | = 0.999 Q. 

The matrix elements of the weak current between the B s meson and the final D s meson 
entering in the factorized nonleptonic decay amplitude (|55|) are parametrized by the set of the 
decay form factors. Using the form factors obtained in Sees. HV ] IVIl rVIIII we get predictions 
for the branching ratios of the nonleptonic B s decays to ground state and excited D s mesons 
and present them in Tables IVHH IIXI in comparison with other calculations and available 
experimental data. We can roughly estimate the error of our calculations within the adopted 
factorization approach to be about 20%. It originates from both theoretical uncertainties in 
the form factors, effective Wilson coefficients and experimental uncertainties in the values 
of the CKM matrix elements (which are dominant), decay constants and meson masses. 

In Table IVTHl we give predictions for the branching ratios of the two-body nonleptonic B s 
decays to the ground state meson and light (ir, p, Ky) or heavy Dy^ meson. We com- 
pare our results with predictions of the QCD sum rules |l8|, relativistic constituent quark 
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TABLE VIII: Comparison of various predictions for the branching fractions of the nonleptonic B s 
decays to ground state D s mesons with experiment (in 10~ 3 ). 





Decay this 


paper 


[18] 


\m 


[20] 


[41] 


[42] [43] 


Experiment [1] 


B s 


-> D-7T+ 


3.5 


5 


9 7+0.2 
z - ' -0.3 


■•"'-0.6 




1.42 ±0.57 1.96±J;|f 


3.2 ±0.4 


B s 


-> D-p+ 


9.4 


13 


6.4±i;? 


4 2+ 1 - 7 




4 7+2-9 
' -2.3 


7.4 ±1.7 


B s 




2.7 


2 


q 1+0.3 
°--0.2 






2.11 ±0.73 1.89±J;|§ 


2.1 ±0.6 


Bs 


-> ^rp + 


8.7 


13 


9.0il;| 






r 90+3^34 


10.3 ± 2.6 


B s 




0.28 


0.4 


o.2i±g : gS 


o.i3l8:8l 




0.103 ±0.051 0.170i8:8e6 




B s - 


-»• z?-k* + 


0.47 


0.6 


n qo+0.05 


0.28l8i 




0.050 ± 0.022 0.28ll8;io9 




B s - 




0.21 


0.2 


n 24+ 002 

u - z ^-0.02 






0.159 ±0.067 0.164i8:8e4 




B s - 


-»■ D*-K* + 


0.48 


0.6 


U - OD -0.07 






0.163 ± 0.086 0.322l8:i2l 




B s 


-> £>7^ 


11 


10 


8-3±i;8 


oc: + 14 
J0 -12 


16.5 


2.17 ±0.82 


5.3 ±0.9 


B s - 




10 


8 








2.62 ±0.93 




B s - 




6.1 


4 


7-0±l;I 






2.54 ±0.57 




B S ^D S 




r 16.1 


12 


is.4±?:8 




24.0 


5.16 ± 1.10 


12.4 ± 2.1 


B s - 


-»• D*-D*+ 


25 


16 


24±| 




31.8 


27.7 ±7.6 


18.8 ±3.4 


B s ^ 




52.1 


38 


47.7±tl 




72.3 


35.0 ±7.8 


45 ± 14 



models [19|, |41|, the light cone [20J and three-point QCD sum rules [421 ] . the perturbative 
QCD approach (43[. Available experimental data [l| are also given. We find reasonable 
agreement between our results, QCD sum rules [l8[ and quark model [Hj predictions and 



experimental data. Results of quark model calculation [41| are slightly larger, while those of 



three-point QCD sum rules [42| and perturbative QCD [43] are slightly smaller. However, 
experimental and theoretical uncertainties are still too large to make possible the discrimi- 
nation between theoretical approaches. 

In Table ILXl we present our predictions for the two-body nonleptonic B s decays to orbitally 
and radially excited D s meson and light or heavy D s meson. They are compared with the 
results of the light cone sum rules [20|, which are available only for decays involving the 
scalar D* s q meson. In general, central values of our predictions for the decays B s — > D* s q M + 
(where M is a light meson) are slightly larger, but both results are compatible within errors. 
On the contrary, for decays B s — > DIqD^* +S) our results are significantly lower, especially 
for the B s —> D*qD+ decay. The same pattern of our predictions and the light cone sum 
rules results [2(| holds also for the B s decays to ground state mesons (see Table IVHT1) . From 
Table IIXI we see that some of the nonleptonic B s decays to the excited D s mesons have 
branching fractions comparable with the ones for the decays to the ground state D s mesons, 
given in Table IVIII1 

Ver y re cently the LHCb Collaboration announced the first observation of the B s — > D s i7r 
decay [44J. Only the relative branching fraction of this decay was measured. However, this 
observation indicates that we can expect the measurement of the nonleptonic B s decays to 
excited Z)„ mesons in near future. 
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TABLE IX: Branching fractions of the nonleptonic B s decays to orbitally and radially excited D, 



mesons (in 10 3 ). 


Decay 


this paper 


[20] 


B s D*qTt+ 


0.9 


0.52t°f 


B s D* sQ p+ 


2.2 


i q+0.6 
—0 5 


B s -> D S VK+ 


0.07 




5 S -> z? s v^*+ 


0.12 


0.08±°g| 


5 S ->• D* s0 Df 


1.1 


13l 5 7 


B s -> D S V-DJ+ 


2.3 


6.0±I1 




0.29 






0.83 






0.021 






0.044 




B s -»• 


0.54 




b s ^1il>:+ 


1.5 




5 S -> DJ l7 r+ 


1.9 




5 S D7 lP + 


4.9 




5 S -> D^K + 


0.14 




B s £>7 1 1T* + 


0.26 




5 S -> £££>+ 


3.0 






5.9 






1.6 






4.2 




B s -> 


0.12 






0.22 






1.4 






4.5 




5 S -> D s (2S)-ir+ 


0.7 




B s -> D s (25)-p+ 


1.7 




B s -»• £>J(25) - 7r + 


0.8 




B a DJ(25)-p + 


2.2 




B s -> D s (25)-K+ 


0.05 




5 S -> D s (25)-i^*+ 


0.08 




5 S -> Z>*(2S')-.K' + 


0.06 




5 S -> L>*(25')-^*+ 


0.12 




B a ->D a (2S)-D+ 


1.0 




B s ^ D S (2S)-D* S + 


0.7 




B S ^D*(2S)-Dt 


0.7 




B s -> D*(2S)-D*+ 


1.7 
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XI. CONCLUSIONS 

The weak form factors of the B s decays to ground state D s mesons, as well as to first 
orbital and radial excitations of D s mesons were calculated in the framework of the rela- 
tivistic quark model based on the quasipotential approach. The heavy quark expansion was 
applied for the calculations of the form factors of the weak B s decays to Dy and (2S) 
mesons. The obtained form factors satisfy all model independent constraints imposed by 
heavy quark symmetry and HQET. The leading and subleading Isgur-Wise functions were 

expressed through the overlap integrals of the meson wave functions. The form factors of 

(*) 

weak B s decays to the orbitally excited D s j mesons were calculated, by using previously 
developed methods 0]. All relativistic effects, including contributions of the intermediate 
negative-energy states and transformations of the wave functions to the moving reference 
frame were consistently taken into account. For the numerical evaluations the relativistic 
wave functions of B s and D s mesons, obtained as the solutions of quasipotential equation 
(JTJ in Ref. [Io[, were used. As a result the weak decay form factors were determined in the 
whole accessible kinematical range without applying any additional parametrizations and 
extrapolations. This significantly reduces theoretical uncertainties of the results. 

Using these form factors we considered various semileptonic B s decays governed by the 
b — > c weak transition. The obtained results were compared with previous calculations based 
on constituent quark models, light cone sum rules and QCD sum rules. The following total 
semileptonic B s branching ratios were found: 

1. for decays to ground state Z?W mesons Br(B s — > D^eu) = (7.4±0.7)% and Br(B s — > 
D&tv) = (1.92 ± 0.15)%; 

2. for decays to radially excited D^{2S) mesons Br(B s ->- D^(23)ev) = (0.65 ±0.06)% 
and Br(B s D^{2S)ru) = (0.026 ± 0.003)%; 

3. for decays to orbitally excited D^*J mesons Br(B s — > D^ev) = (2.1 ± 0.2)% and 
Br(B s -)• D { *)rv) = (0.11 ± 0.01)%. 

We see that these branching fractions significantly decrease with excitation. Therefore, 
we can conclude that considered decays give the dominant contribution to the total 
semileptonic branching fraction Br(B s — > D s ev + anything). Summing up these contri- 
butions we get the value (10.2 ± 1.0)%, which agrees well with the experimental value 
Br(B s -> D s eu + anything) Exp . = (7.9 ±2.4)% Q. Note that our predictions for the branch- 
ing ratios of semileptonic decays to orbitally excited states B s — > D s \\iv and B s — > D s2 pv are 
in reasonable agreement with recent data from the DO [33[ and LHCb [HJ Collaborations. 

The tree-dominated two-body nonleptonic B s decays to the ground state or excited D s 
meson and the light or charmed meson were calculated in the framework of the factorization 
approximation. This allowed us to express the decay matrix elements through the products 
of the weak form factors and decay constants. The obtained results were compared with 
previous calculations and experimental data, which are mostly available for the decays in- 
volving ground state D s mesons. Good agreement of our predictions and data was found. 
Detailed predictions for decays involving orbitally D^J and radially D^*\2S) excited mesons 
were obtained. Some of such decays have branching fractions comparable with the ones for 
decays to ground state D s mesons. The following decay channels were found to be the 
most promising: (1) decays to excited D s and light mesons B s — > D^p + , B s — > D* S 2P + , 
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B s -> L>: Q -p+, £?, -> D*(2£)-p+ 5 S -> L>7 l7 r+, £? s -> D s (2S)-p+, S s -> (2) de- 

cays to excited and ground state D s mesons B s — > D~i_D* + , B s — > Z)*^.D* + , £> s — » D~ X D+ . 
Therefore there are good reasons to expect that these decays will be measured in the near 
future. This expectation is confirmed by the very recent observation of the B s — > D s in 
decay by the LHCb Collaboration jUj]. 
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Appendix A: HQET expressions for the weak form factors of the B s decays to 

ground state D s mesons 



In HQET the weak form factors of the B s decays to ground state D s mesons up to 1/ rriQ 
order are expressed as follows |l5| 

(Al) 



K = £ + {e c + e h ) [2 X i - 4(iu - 1) X 2 + 12* 3 ] , 



e c -e b ) 2£ 3 - A£ 



(A2) 



£ + e c 2 X i - 4 X 3 + M 



2 X i - 4(w - 1) X 2 + 12 X 3 + M~ 26 



h Al = £ + e 



w 1 - 

2xi - 4X3 + — — rA£ 
w + 1 



2xi - 4{w - l)x 2 + 12X3 + ^4 ~ 2^ 



w + 1 



h>A 2 — £c 
hA 3 



4X2 



w + 



W 1 - 

2xi - 4X2 - 4x3 + —^rM 



w + 1 w+1 
+e b [2 X i - i{w - 1)X2 + 12X3 + M- 26 
where Eq = 1/(2ttt.q) and A = M — uiq. 
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(A3) 



(A4) 
(A5) 

(A6) 



Appendix B: Relations between two popular sets of form factors 



M Ba + M Dau (Ml+Ml-q 2 



M 



2M B M Da 



u,-M D fM k + M k ^\ 



2^M Bs M Da 



2Mb Md. 



Hq 2 



2^M Ba M Ds 



(M Ba + M Da f - q\ h (Ml + Ml - q^ 



Mb. + M Da 



2M B Mn 



26 



(M Bs -M Ds f-q\ f M 



2M Ba M Da 



M Bs +M D * fM 2 +M 2 -q 2 



V(q 2 ) = ' h v \ 

2yfM Bs M D . V 

(M Bs + M Dt ) 2 -q 2 



2M Ba M D% J 

M 2 +M 2 -q 2 ^ 



2JM B M D JM Bs + M DsJf ) 



-h 



2M B M D% 



A,{q 2 ) = 



M Bs + M D% 
2^M Ba M D% 



( Ml+M^-f \ Mp* (M Bs +Ml t -q* 
h M — ^J-JJ-. — + T7T h M - 



\ 2M B M D% 



M Ba "\ 2M B M Dt 



Mq)= - " H 2M B M m , 



4M D|V 


lM B M D * 


Ml 


-M 2 * \ 

^ s 


4M D .y 


/M Bs M D * 


+ 


^ \ 


AM D * X 


lM Ba M D% [ 



hA A — 9aZT~/i7Z — + irr h M " 



\ 2M B M D - t 



M Ba "\ 2M B M Dt 



h A 3 [ 



fM 2 +M 2 -q 2 \ M D . fM 2 +M 2 -q 



\ 2M Ba M D , 



hA 2 \ 



M Ba "\ 2M Ba M D% 



(B2) 
(B3) 
(B4) 
(B5) 



, (B6) 



(*) 

Appendix C: Helicity components of the hadronic tensor for the B s — > D y s 'lu decays 

(a) B s — > D s transition 



H± = 0, 

A V2 



H 



i 



(M 2 Bs -M 2 Da )f (q 2 ). 



(b) B s — > D* transition 

A 1 / 2 



H (q 2 ) = 



M Ba + M D 
1 



V(q 2 ) T 



(M Ba + M D * 
A 1 / 2 



(M Ba + M Dt ){M 2 Ba - M 2 - q^A^q 2 ) 



A 



M Ba + M D * 



A 2 (q 2 ) 



\ 1/2 

H, = —A (q 2 ). 



(CI) 



(C2) 
, (C3) 
(C4) 



Here the subscripts ±,0, t denote transverse, longitudinal and time helicity components, 
respectively. 
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Appendix D: HQET expressions for the weak form factors of the B s decays to 

radially excited D s mesons 

In HQET the structure of the weak decay form factors for B s decays to radially excited 
D s [(n + l)^] mesons up to 1/rriQ order is the following 27 



K 

fi- 
tly 



b A 3 



e (n) + e c [2xi - 4(w - l)x 2 + 12* 3 ] + e bXb , 



w — 1 



26 



+e b 

£ {n) + e 
+£b jx& + 
e c |4X2 - 

i {n) + e 



u> — 1 



2xi - 4xs 
w — 1 



w — 1 

w — 1 



iu + 1 
2 



A 



w + 1 

AW - A 

w — 1 



w — 1 



W + 1 

2xi - 4x 2 - 4x3 + 



it; — 1 
u; — 1 



ty + 1 



A (n) + 



A(«) _ 

1(7 — 1 



IU + 1 



6 



(Dl) 
(D2) 



(D3) 



(D4) 
(D5) 



(D6) 



where A(AW) = M(M^) — mg is the difference between the heavy ground state (radially 
excited) meson and heavy quark masses in the limit tuq — > oo. 



Appendix E: Helicity components of the hadronic tensor for the B s 



(a) B s —y D* transition 

H± 
H 

H t 

(b) B — > D s \ transition 

H± = (M Ba + M Dsl )h Vl (q 2 ) ± 
1 



0, 

A 1/2 

7¥ 
1 

7¥ 



[(M 2 Ba -M 2 Ds0 )r + (q 2 ) + q 2 r4q 2 )]. 



X 1 ' 2 



-h A (q 2 ), 



Hn 



2M DslV ^ 



M Bs + M Dal 
T (M Bfl + M Dsl )(Mi s - M 2 si - g 2 )^(g 2 ) 



D^*jlu decays 



(El) 
(E2) 

(E3) 



(E4) 
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A 



H, 



2M Bs 

A l/2 



[hv 2 (q 2 ) + h V3 (q^, 

(M Bs + M Dal )h Vl (q 2 ) + ^ [Mg 2 ) + M<? 2 )] 



(E5) 



2M 



2M 



■[M<? 2 )-M<? 2 )] • 



(E6) 



(c) B — > D' sl transition 

Hi are obtained from expressions (1E4I) - (1E6I) by the replacement of form factors hi(q 2 ) by 
gi{q 2 ) and the final meson mass Md s1 by . 

(d) -B — >■ D* 2 transition 



H 



A 1 / 2 



H, 



2V2M Ba M Da2 
A 1 / 2 

2^M Ss M 2 s2 v^, 

+-^itA 2 (q 2 ) + t A3 (q 2 )} 

[2 A 



(M Bs +M Ds2 )t Al (g 2 )± 



A V2 



-t v {q 2 



M Bs + M Ds2 
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(E7) 



(E8) 



V3 4M Bs M 2 s2V ^ 
+ 2^ [ ^ (g2) -^ (g2)] 
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2M 
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